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Electron-electron interactions strongly affect the behavior of low-dimensional systems. In 
one dimension (ID), arbitrarily weak interactions qualitatively alter the ground state pro- 
ducing a Luttinger liquid (LL) 1 which has now been observed in a number of experimental 
systems 2-6 . Interactions are even more important at low carrier density, and in the limit 
when the long-ranged Coulomb potential is the dominant energy scale, the electron liquid is 
expected to become a periodically ordered solid known as the Wigner crystal 7 . In ID, the 
Wigner crystal has been predicted to exhibit novel spin and magnetic properties not present 
in an ordinary LL 8-12 . However, despite recent progress in coupled quantum wires 1314 , un- 
ambiguous experimental demonstration of this state has not been possible due to the role 
of disorder. Here, we demonstrate using low-temperature single-electron transport spec- 
troscopy that a hole gas in low-disorder carbon nanotubes with a band gap is a realization of 
the ID Wigner crystal. Our observation can lead to unprecedented control over the behavior 
of the spatially separated system of carriers, and could be used to realize solid state quantum 

computing with long coherence times. 

Carbon nanotubes are high mobility quantum wires that may enable the study of the intrinsic 

properties of the ID electron gas without interference from disorder. Single quantum-dot transport 
experiments have been performed, demonstrating Coulomb blockade 1516 and Kondo physics 17 , 
down to the few-electron-hole regime 18-20 . These experiments have generally been interpreted 
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using the orthodox model and four-fold shell- filling. However, deviations have been noted at 
low density 18 20 , suggesting the importance of electron-electron interactions. In ID, the ratio of 
inter-electron Coulomb energy to kinetic energy goes as l/(na B ), where n is the carrier density, 
a B = eh 2 /me 2 is the Bohr radius, e is the dielectric constant, e is the electric charge, h is Planck's 
constant, and m is the electron effective mass. On the basis of this order-of-magnitude estimate 
(see e.g. ref. 9), the Coulomb energy should dominate and the system should begin to cross over 
to a strongly interacting regime at n ~ as' 1 = me 2 /eh 2 . While most previous studies in the 
few-electron regime have been done on nanotubes with small bandgaps (<100 meV), observing 
the Wigner crystal state is more favorable in nanotubes with larger gaps, where the effective mass 
is larger. 

We report axial magnetic field measurements on low-disorder carbon nanotubes with larger 
bandgaps than previous studies. We observe for the first time three distinct regimes as a function of 
carrier density and magnetic field: (I) a completely spin and isospin polarized state, (II) an isospin 
polarized, spin antiferromagnetic state, and (III) an unpolarized state with a four-fold addition en- 
ergy period. The transitions among these regimes can be quantitatively and intuitively explained 
using a Wigner crystal picture based on a gapped LL model 21 with the carriers represented by spa- 
tially localized solitons. The value of the soliton width we obtain from our analysis, corresponding 
to the spatial extent of the carrier wavefunctions, agrees with theoretical estimates 2122 . Surpris- 
ingly, we observe a sudden quenching of the Kondo effect for odd hole states as the magnetic field 
is tuned through the transition from the antiferromagnetically ordered regime to the four-fold fill- 
ing regime. This can be understood based on the formation of large spin states with spin S > 3/2 
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resulting from the interplay between the spin and isospin degrees of freedom in the Wigner lattice. 

Nanotubes with a bandgap are known to be more susceptible to disorder than metallic nanotubes 23 . 
Our fabrication procedure, based on Cao et al. 20 , is designed to reduce disorder by growing nan- 
otubes in a suspended geometry to eliminate perturbations from the substrate and directly over 
metal contacts in order to avoid disorder due to chemical processing. Fig la shows a scanning 
electron microscope image of a representative device, while Fig lb shows a schematic diagram of 
our devices, consisting of a nanotube with attached source, drain and gate electrodes. 

At low temperatures, we typically observe Coulomb peaks in the conductance G versus the 
gate voltage V g , when each additional hole is added to the nanotube dot. Fig lc shows such G 
versus V g data taken at temperature T = 1.4 K from a device Dl with length 500 nm. V g >1 
V depletes the nanotube, placing the Fermi level in the bandgap. The regularity of the Coulomb 
peaks indicates that we are studying a single quantum dot that is essentially free from disorder 
and pinning. We also perform non-linear transport spectroscopy; Figure 2a shows a color-scale 
plot of dl/dV versus source-drain bias V and V g for Dl. The data exhibits Coulomb blockade 
diamonds, shown by the dashed lines, and analysis of the data yields a charging energy U ~10 
meV, consistent with that expected for a 500 nm-long nanotube dot. We also measure the scale 
factor a between energy and gate voltage and find for Dl a « 14. Using the value for a and a 
measurement of G versus V g over the bandgap, we infer the bandgap 2 A, which for Dl was found 
to be ~220 meV. 

Each added hole has both a physical spin and an isospin corresponding to the sense of the 
orbital motion of the carriers around the tube waist 19 . In a magnetic field B, the ground state 
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energy increase from each hole shifts by the sum of contributions E orb from its orbital magnetic 
moment 19 - 20 and the Zeeman energy E z from the hole spin 24 . These energies are given by E orb = 
±fi orb B = ±rev F B/2 « 0.41meVB[T]r[nm] 19 ' 20 , and E z = ±g^ B B/2 « ±0.0585[T] meV, 
where fi B is the Bohr magneton, r is the tube radius, e the electric charge, v F the Fermi velocity, 
g « 2 is the electron ^-factor. As the position of iVth Coulomb peak in V g is a factor a times the 
energy difference between the N and iV + 1 hole ground state, we determine the spin and isospin 
for each added hole by applying B parallel to the tube axis (Fig. 2b) and studying the Coulomb 
peaks' shifts versus B. Our results are plotted in Fig. 2c, which shows a color-scale plot of G vs. 
B and V g from Dl. 

Three distinct regimes are evident. The peaks up to number ~ 12 are nearly parallel from 
near B = T up to B = 8 T, shifting to larger V g for larger B (positive slope). This region 
is referred to as region I. At more negative V g , the slopes alternate between two different positive 
values (region II). Finally, a region where both positive and negative slopes are observed, as well as 
approximately horizontal features occurring quasi-periodically in B (region III). The abrupt change 
of the slopes of the Coulomb peaks, particularly evident at the region II-III boundary indicates a 
change in the ground state wavefunction configuration. Curves are shown superposed on the graph 
approximately delineating the three regions corresponding to the results of a theoretical model fit 
that will be discussed below. 

Figure 3a shows the measured slopes at B = 8 T in region I and II versus peak number. The 
peaks in region I have nearly the same slope, with d(E z + E orb ) jdB ~ 0.4 meV/T, while peaks in 
region II show a pronounced alternation in addition energy, with an amplitude AE S ~ 0.14 meV/T, 



4 



and generally falling below the data in region I. This indicates that in region I, all the holes enter 
the dot with the same spin and isospin. We thus find fi orb ~ 0.33 meV/T, yielding an estimated tube 
radius of ~0.8 nm. In region II, the holes enter with alternating spins, but with the same isospin. 
From AE S we estimate the electron g-factor to be «2.4, in good agreement with the theoretical 
value. Similar behavior has been found in all devices with a significant gap that showed regular 
Coulomb oscillations, attesting to the generality of this behavior in low-disorder samples. 

The spin polarization observed in region I is surprising. In a shell-filling picture, the subbands 
corresponding to different isospin are split by the field 19 ' 20 , and holes first begin to fill the lowest 
energy subband. This could account for the consistent isospin of the added holes. Yet, we expect 
because of the spin degenerate hole wavefunctions, they should enter with alternating spins except 
at sufficiently high magnetic field B c where the Zeeman energy would exceed the mean level 
spacing AE. Taking a parabolic potential as in ref. 18, we find AE ~ 1.5 meV, giving B c = 13 
T, precluding entirely the existence of region I over the range of magnetic fields studied. Using a 
hardwall potential, as done in ref. 20, but with the additional inclusion of spin yields significantly 
fewer electrons in region I than we observe and also predicts for region II level crossings that would 
produce kinks in the evolution of the Coulomb peaks and singlet-triplet Kondo-valley conductance 
enhancements 25 , which we do not observe. Thus, we cannot readily account for the observed 
behavior using simple shell-filling pictures, consistent with the strong interactions expected at 
these low densities. 

This behavior, however, can be naturally accounted for by assuming that the holes form a 
ID Wigner crystal, in which the carriers become periodically ordered. Considering spin only, the 
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inter-carrier exchange coupling J 9 is expected to be antiferromagnetic, as opposite spin carriers 
can both occupy their lowest energy molecular orbital. In nanotubes, the additional degree of 
freedom provided by the isospin results in modifications to the usual picture of exchange based on 
spin alone: the minimum energy state occurs when nearby carriers have different spin or isospin 
quantum numbers. Otherwise, the nearby carriers pay an energy cost J. When n increases and 
the carrier wavefunctions overlap more strongly, a four fold addition energy period is expected; 
after exhausting all four spin and isospin combinations, adding the next carrier requires additional 
exchange energy. 

When B > the state of the Wigner crystal is a competition between magnetic and exchange 
energy. Since J is predicted to become exponentially small at low n n ' 26 , we expect at low n the 
total energy will be minimized by having the holes both spin and isospin polarized. At intermediate 
n, since » E z , we expect a transition to an antiferromagnetically spin ordered, isospin 
polarized chain when gfXsB = 2 J. For larger n, where non-nearest neighbor exchange becomes 
important, we expect to observe a four-fold addition energy pattern as discussed above. This 
picture thus qualitatively accounts for the existence and behavior of the three distinct regions as 
well as the transitions among them as n increases. 

To quantitatively analyze our data, we utilize the theory of Levitov and Tsvelik 21 . The carri- 
ers in the ID Wigner crystal within the nanotube are treated in terms of a linear combination of a 
charged bosonic field and three neutral flavor bosonic fields, chosen to yield their given spin and 
isospin. The fields' behavior is governed by a kinetic energy cost for localization that competes 
with a gap-dependent potential energy gain. This yields spatially-localized soliton solutions for 
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the fields, with an optimal width w for the flavor solitons. Theoretical estimates yield a range of 
w ~ 5 — 20 nm for sample Dl 21 - 22 . Figure 3b shows the classical spatial variation of the four fields 
necessary to produce a chain with the four different spin and isospin combinations. 

We use the theory of Levitov and Tsvelik in a classical approximation (details are provided 
in the Supplementary Discussion) to compute the ground state energy for a spin polarized (I), 
antiferromagnetically ordered (II), and four-fold period (III) ground state using the spin and isospin 
configurations shown in the insets to Figs. 3a and 3b. As n increases, the flavor solitons overlap 
more strongly, increasing J. Thus, we expect as n increases the ground state for B > will make 
transitions among the three states. However due to e.g. the non-uniform potential along the tube, 
we expect that the transitions between these states should happen via a sequence of individual spin 
or isospin flip events with a similar energy per carrier as calculated by our method. Using toasa 
fitting parameter, we fit the observed transition between regions II and III to the calculated field 
where the two states are degenerate. From the data in Fig. 2c we find w ~ 9 nm. This fit is 
superposed on the same data as the yellow curve and follows the observed boundary closely. The 
fitted value agrees well with the range of theoretically expected values. 

Using the fitted value for w = 9 nm, the boundary between regimes I and II is computed 
directly, plotted in Fig. 2c as the red curve. This curve indicates the expected number of spin- 
polarized holes added at B = 8 T is ~9 in agreement with the observed number ~12, although 
this boundary is not as distinct as that between II and III, which we attribute to thermal fluctuations 
that should blur the boundary by ~ k B T/g^ B ~ IT. Nevertheless, given that the curve is plotted 
without free parameters, the agreement is satisfactory. 
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Finally, analysis of the transition between regions II and III reveals a unique aspect of Wigner 
crystallization. This is shown as a color-scale plot of the conductance, G versus B and V g in Fig. 
4a. Regions II and III are shown, divided by a black dotted line. In region II, G between adjacent 
Coulomb valleys alternates, with larger G in the odd-hole valleys, as shown in Fig 4b. This is a 
well-known signature of the Kondo effect 27 , which results from the screening of a local spin on the 
nanotube with the conduction electrons in the electrodes to form a spin singlet state. This state has 
characteristic binding energy ~ k B T K , where T K is the Kondo temperature 27 . The conductance 
alternation between valleys, and slopes of the Coulomb peak motion discussed previously indicates 
that the spin of the nanotube changing between S = S in the even valleys and S = S + 1/2 in 
the odd valleys, where So is a constant 24 . 

Remarkably, when B is tuned from region II to region III in an odd-hole valley (from A to B 
Fig. 4a), G drops by AG ~ 1.5 fiS. In contrast, in the even-hole valley no such drop is observed 
(Fig. 4a points C and D). Fig. 4c shows G versus B taken along the red and yellow vertical dashed 
lines in adjacent Coulomb valleys. 

A shell-filling model could possibly account for the ridges by level crossings and an enhance- 
ment of T K because of the orbital degeneracy 28 . However, this picture cannot readily account for 
the observed conductance drop AG since an odd hole state with total spin S — S + 1/2 should 
remain in the same spin state even when the up and down isospin bands first align at the Fermi 
level, as illustrated in Fig. 4d. This picture therefore predicts AG = 0, contradicting the observed 
dip. 

In a Wigner crystal, the transition between region II and III occurs when B is lowered so 
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that a hole flips its isospin. This hole also flips its physical spin to minimize its magnetic energy 
without an exchange energy cost. The total spin state becomes S — S + 3/2, as shown schemati- 
cally in Fig 4e. Since T K generally becomes lower as S increases 25 , G should drop between A and 
B, as observed. Thus the observation of AG > in the odd-hole valley provides a clear exper- 
imental signature of the ID Wigner crystal state. Taken together with the quantitative agreement 
to the theory of Levitov and Tsvelik, our data provides convincing evidence that carriers in carbon 
nanotubes at low densities form a ID Wigner crystal. 

The realization of this long predicted state can now be used to test theories of interacting 
electrons in ID in the clean limit 9-12 - 22 . For instance, at B = the equilibrium state is determined 
by a competition between the thermal energy k B T and J. If k B T » J, the spins and isospins 
can flip freely because of thermal fluctuations. This spin-incoherent regime is predicted to exhibit 
different behavior from an ordinary quantum wire, for example reduced conductance, which may 
be related to the 0.7 structure observed in quantum point contacts 9 . 

Moreover, we note that at low density, the experiment achieves a carrier separation of ~ 100 
nm which gives experimental access to control of individual exchange couplings, say using local 
gates. One could then utilize the many-body Wigner crystal as a chain of quantum bits towards 
realizing spin-based quantum computing in carbon, where intrinsic spin lifetimes are expected to 
be longer than in conventional semiconducting materials. 
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Supplementary Discussion 



Calculation of exchange energy 



The Levitov and Tsvelik theory 21 treats the carriers in a nanotube as a gapped Luttinger 



liquid (LL). The interactions in the gapped LL are characterized by a charge stiffness parameter 
K which is related to the inter-carrier interaction potential and can be usefully approximated by 
a constant ~ 20 — 40 21 ' 22 . The interactions renormalize the non-interacting electron bandgap, 
making it ~ K 1 / 2 times larger. 

In a nanotube LL, the two spin and two isospin states are described by one charged and three 
neutral flavor bosonic fields. When a spectral gap is introduced into the LL, the minimum energy 
field conguration becomes a lattice of sine-Gordon solitons, concentrating charge, spin, and isospin 
into localized packets to form the ID Wigner crystal. The characteristic width w c of the charge 
solitons is related to the measured gap A by w c ~ hvj A 29 , where v c = K x l 2 v F is the charge 
mode velocity 21 . The flavor soliton width is w ~ K~ l l 2 w 2 ^ 22 . 

We compute the exchange energy in the Wigner crystal starting from the bosonized La- 
grangian for the gapped Luttinger liquid with the four bosonic fields, which is the sum of two 
terms V + V\. The competition between V and V\ leads to soliton solutions for these fields. 
The solutions 21 , together with the potential, enable us to compute the energy treating the fields 
classically. 

The first term V arises from the kinetic energy of localization 1 ' 21 



where O = 4>c+ is the charged field 01,2,3 = 4> s +-,4>s--,4>c- are the three flavor fields, q is a 




a=l 



(1) 
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wavevector, a; is a position coordinate, v F is the Fermi velocity, and k is Plank's constant. The 
second term V\ is an effective potential that arises from the gap, given by 21 

V 1 = -4A J dx[cos((f) c+ ) cos(0 c _) cos(0 s +) cos(0 s+ ) + sin(0 c+ ) sin(0 c „) sin(0 s+ ) sin(0 s+ )], 

(2) 

Where A ~ %vf/{^w 2 ) is a (renormalized) coupling constant that is related to the physical flavor 
soli ton width w. 

A single soli ton centered at the origin for the charge field is approximated by 21 

f (x) = j [e^6(-x) + (2 - e~^) 9(x)] , (3) 

where 9(x) is the unit step function. As a composite soliton of charge and flavor modes is the 
lowest energy state 21 , the flavor fields then vary as 

/<(*) = (3~ [eZe(- x ) + (2 - e -S) 0(x)} , (4) 

Where j3={ 1,1,1) for a spin up, isospin up carrier, f3=( 1,-1,-1) for a spin up, isospin down carrier, 
f3={- 1,-1,1) for a spin down, isospin up carrier, and j3={- 1,1,-1) for a spin down, isospin down 
carrier. A soliton lattice corresponding to a chain of N carriers with density n is given by 

N 

M x ) = J2fo( x ~j/ n ) ( 5 ) 

N 

<t>i{x) = ^2Pijfi(x - j/n), (6) 

where are the numbers (3i for the jth carrier. The exchange energy is computed by substituting 
the expressions for <fi (x) and <fii(x) into the the potential V + V 1 and computing the energy nu- 
merically. The magnetic energy is computed by multiplying the net spin by g/j, B B/2 and the net 
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isospin by (J> or bB. The sum of the exchange and magnetic energy for a spin and isospin polarized 
electron or hole lattice, a spin- alternating, isospin polarized lattice, and a four-fold period lattice is 
used to determine which of these three states is the ground state. The result for the boundary be- 
tween the four-fold and spin- alternating region is then converted into a mathematical interpolation 
function to fit to the data. 

We note that this model is expected to be valid as long as the soliton width significantly 
exceeds the tube diameter 21 , as it is expected to theoretically, and as it does in our fitted results. 
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Figure 1 Experimental geometry and characteristic transport data, a, Scanning elec- 
tron microscope image of a representative suspended nanotube device. Samples are 
fabricated by chemical vapor deposition growth of single-walled carbon nanotubes, from 
lithographically defined Fe-Mo catalyst islands 30 across predefined Pt electrodes 20 , by 
flowing a gas mixture of methane (0.5 SLM) and hydrogen (0.7 SLM) at 800 °C for 5 
min. Devices are first studied at room temperature in an inert environment and the gate 
voltage and high bias (up to 1 .5 V) electrical transport characteristics are recorded. De- 
vices determined to have a bandgap and that show negative differential conductance at 
high bias with a maximum current of ~10/L /iA (where L is the known device length in 
/iim), corresponding to individual, suspended, single-wall nanotubes 31 , are then selected 
for low-temperature measurement, b, Schematic of device geometry showing nanotube 
with attached source, drain, and doped Si gate electrode ~600 nm below the nanotube. 
c, Conductance versus gate voltage for a 500 nm long device at temperature T=1 .4 K. 

Figure 2 Transport spectroscopy and magnetic field evolution of Coulomb peaks, a, 
Color plot of dl/dV versus gate voltage V g and source-drain bias V. b, Magnetic field 
B orientation relative to nanotube. c, Color plot of G versus V g and magnetic field B. 
The three different regions and curves, obtained from a model fit, delineating them are 
described in the text. We expect similar behavior in the electron-doped regime due to 
electron-hole symmetry 18 , but the Pt electrodes generally produce Schottky barriers for 
electrons that make it difficult to obtain very clear data for the electron-doped regime in 
larger gap devices. Nevertheless, such data should be possible to obtain in principle, 
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perhaps using a different contact metal. 

Figure 3 Energy shift of Coulomb peaks with B and schematic diagram of solitons cor- 
responding to the four combination of spin and isospin. a, Energy shift per tesla dE/dB 
of Coulomb peaks measured at B=8 T in region I and II. dE/dB is roughly constant in 
region I and shows alternation in region II. The value of the tube radius we find from anal- 
ysis of this data is consistent with the range of radii typically reported for methane-based 
nanotube growth over the catalyst we employ 30 , although knowledge of the tube radius 
is not essential to our analysis and interpretation. Note that the radius does not have 
to be related to the measured gap in general, since gaps can be opened due to radius- 
insensitive mechanisms like axial strains and twists 19 and are expected to be enhanced 
due to interactions 2122 . Note also that the results from region I differ from those of Tans et 
al. 32 , in which repeated addition of the same spin electron to a metallic nanotube was ob- 
served for B <~ 1 T, followed by a complex spin-filling pattern at larger B. Here, the spin 
polarization is observed in a gapped nanotube only at low densities and over the entire 
measured range of B up to 8 T. Inset: schematic diagrams of spin configurations in region 
I and region II. The holes all have isospin up. b, Spatial variation of the four fields corre- 
sponding to charge sum in the two isospins (0 C+ ), and three flavor modes corresponding 
to the charge difference (0 C _), and spin sum (<f> s+ ) and difference (0 S _). The correspond- 
ing hole configuration is shown at the top with the arrows indicating the spin and the color 
indicating the isospin, with blue indicating isospin up and red, hatched isospin down. The 
flavor soliton width w is indicated by the red lines. The solitons are here depicted as non- 
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overlapping. When there are ~ L/w ~ 50 holes in D1 the solitons should merge together. 
This corresponds to a nearly constant charge and flavor density. In this limit, we expect 
that the device should therefore cross over to a metallic hole-liquid like regime, where the 
behavior should be well-described by a conventional shell-filling picture, consistent with 
the well-defined four-fold periodicity of the Coulomb peaks in this regime. 

Figure 4 Kondo effect and flavor configurations, a, Color scale plot of G versus V g and 
B. The boundary between regions II and III is shown by the black dashed line. Points 
A and B are in an odd-hole valley, while points C and D are in an even-hole valley, b, G 
versus V g showing Coulomb peaks in region II taken at B = 5 T. c, Line traces from part 
a taken along the odd-hole valley (shown in red) and in the even-hole valley (shown in 
dark yellow), d, Shell filling picture of the transition between region II and III. Isospin up is 
shown by blue levels, while isospin down is shown by red levels, e, Wigner crystal picture 
of the transition between region II and III. 
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